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Abstract 

A non-iterative method is presented for the factorization step of sector decompo- 
sition method, which separates infrared divergent part from loop integration. This 
method is based on a classification of asymptotic behavior of polynomials. The prob- 
lem is converted to ones for convex body in Euclidean space. They are solved with 
algorithms developed in computational geometry. A test implementation shows that 
this method produces less number of decomposed sectors than usual iterative sector 
decompositions. 
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1 Introduction 

In the calculation of Feynman amplitude with massless particles, one has to regulate infrared 
divergences (IR), which cancel out among loop corrections and real emission processes. In 
perturbative QCD, this divergences are regulated by D-dimensional method. Divergent part 
is expressed as poles in terms of e = (4 — D)/2. Sector decomposition method developed in 
Refs. jl] 12 |4] are widely used for separating IR divergences. 

As a simple example of separation of infrared divergences, let us consider the following 
one-dimensional integral for e < with a regular function / such that /(0) is a non-zero 
finite value: 

1=1 dxx- l - e f{x) 







1 dx x-^fiO) + [ dx x- 1 -* (f(x) - /(0)) (1) 
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When e = 0, the factor x _1 ~ c produces a logarithmic divergence by the integration around 
x ~ 0. While for e < 0, the divergence is regularized and converted to a pole of e. The 
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coefficient of the pole is /(0), which is the first term of the Taylor expansion of f(x) around 
x = 0. The rest of the integration becomes finite for e->0. This example shows that, when 
the singular part of the integrand is factored out, one can separate divergent part as poles 
in terms of e. 

For the general loop integrations, it will also be possible to separate IR divergence, if 
the singular part is factored out. Although the factorization of a multi-variate polynomial 
is not a trivial problem, it is solved by the method of sector decomposition. This method 
shows that when the integration domain is properly decomposed into sectors and selecting 
good variables in each sector, the divergent part is factored out. This method provides at 
the same time a practical procedure to find such sectors and variables. Then we can obtain 
the coefficients of the poles of e and the finite part as a Laurent series in terms of e. 

Usual sector decomposition method divides the integration domain and find appropriate 
new variables by an iterative way. It was found that a simple iterative method may fall into 
an infinite loop and strategies are proposed to avoid this problem by Refs. [5J E]. Other 
strategies have been proposed and a practical system has been constructed by Refs. [7JIEJIH]- 
Since one has to calculate integration for each sector, a method is preferable when it produces 
less number of sectors. 

We propose another method based on the classification of asymptotic behavior of poly- 
nomials around the origin. The problem is converted to one for convex bodies in Euclidean 
space and solved with algorithms developed in computational geometry. These are deter- 
ministic algorithms without iterations. This talk presents, based on Ref. |10j . our basic ideas 
and show the results of this method. 

2 Sector decomposition 

The procedure of sector decomposition consists of the following steps: 

1. Start from Feynman parameter representation of a loop integration. 

2. Primary sector decomposition: 

The ^-function is integrated and the integration domain is decomposed such that the 
singularities appear in the specific position in each sector. 

3. Factorization of the integrand: 

The sectors obtained in the previous step are decomposed into finer ones. New variables 
are found such that singular parts are factored out. 

4. Separation of poles in terms of e: 

The regular part of the integrand is expanded in terms of e. The coefficients of Laurent 
series in terms of e are obtained. 

5. Integration of coefficients: 

The coefficients are now free from IR divergences. However, they will still include 
multi-dimensional integration. They will be integrated out by analytic or numerical 
methods. 

Feynman parameter representation of a loop integration includes ^-function which defines 
the (N — l)-dimensional hyperplane. The primary sector decomposition integrates over one 
of integration variables and decompose the integration domain into N sectors. An example 
is shown by Fig. [T] for the 3-dimensional case. By the decomposition of the integration 
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(a) 



(b) 



Figure 1: An example of primary sector decomposition in 3-dimensional case. Integration 
plane is divided into 3 sectors and sub-domain abed in (a) is mapped to a square shown by 
(b). 



domain and selecting new variables, sectors become (N — l)-dimensional cubes. The original 
boundary of the integration domain is mapped to a part of the new boundary specified by 
conditions that some of new variables are 0. The resulting expression of the integration for 
l-th sector becomes: 



Jo 



N-l 



Ht v ~ l 



(2) 



where Ti{t) and Ui{t) are polynomials of variables t. Parameter Vj is integer determined 
by the power of propagators, and I is an (N — l)-dimensional vector whose elements are 
all 1. Powers j5 and 7 are functions of e. If Fi{t) or Ui(t) become at the boundary, and 
corresponding power /? or 7 becomes negative integer for e = 0, the integration diverges. 



3 Geometric method 

We want to express polynomial Ti{t) in the form: 

F l {t{z)) = C a z b °{l + H a {z)), H a (0) = 0, z K =l[zf a) \ (3) 

3 

where z is a set of new variable, H a {z) is a polynomial of z, b a is an integer vector and C a 
is a constant. In order to seek such an expression, we consider the following example taken 
from one-loop box integration: 

Fi{t) = -s 23 t 2 t 3 - si2*i - s 4 iit 3 = t 2 t 3 (-S23 - S12T-7 s iT ) • ( 4 ) 

If terms t\jt 2 tz and t\/t2 — > when t 2 H — > 0, the asymptotic behavior of Ti is determined 
by the term —523^2^3- The singular behavior of the integrand around the origin is determined 
by this term. This condition is satisfied by taking new variables z defined by t 2 = z 2 , t 3 = 
z 3 , t\ — z\z 2 z^. With the range of Zj being limited to (0, 1), a sub-domain of the integration 
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domain in i-space is obtained. This example shows that a term of a polynomial in some 
sub-domain will determine asymptotic behavior of the polynomial around the origin. We 
call this term dominant in this sub-domain. 

In order to see dominant terms more closely around the origin, let us change variable to 
Uj = —\ogtj or tj — e~ Vi . Monomial t c — t^t^ 2 •••t^T 1 I becomes e~^ c ' v \ where (c, y) = 
J2 c jVj i s the inner product defined in (N — l)-dimensional Euclidean space. This shows 
that a monomial corresponds one-to-one to a integer vector c in this Euclidean space (let 
us call this power vector). Let Z^ 1 is the set of power vectors corresponding to the terms 
included in a polynomial T\. Polynomial Ti is expressed by: 

Tl{t)= a ct c = /Z a ce-^ y) . (5) 

In order to find the dominant term, let us consider a limit A — > +00 for yj = Xuj — > +00 
with a fixed non-negative real vector u. A term t b = e~^ h ' v ^ is dominant in this limit when 
(b, y) < (c, y) for all c G Z^ 1 . Conversely, let us fix vector b and vary u or y. The term with 
this power vector b is dominant in the sub-domain defined by: 

Af :={y6R^- 1 |( C -6, 2 /)>0,VceZ^}. (6) 

This sub-domain forms a convex polyhedral cone. The polynomial Ti is expressed as the 
following: 

m(v))= E 9 (y £ Af'K ( ^ ) [« b + £ a c e-^~ b ^]. (7) 

bez^i ' cez^i-ib} 

For a vector y e A^ ! , factor e~( c ~ b ' y ) — > for \y\ — > 00, since (c — b,y) > holds for all 
c. This implies that the term t b is dominant in the sub-domain A^ 1 . Thus, asymptotic 
behavior of Ti is classified in terms of these sub-domains {A^ ! }. 




Figure 2: An example of geometric sector decomposition 



When this method is applied also to Hi, the singular parts are factored out from the 
integration. Let us return to our example Eq. ^ of one-loop box. In this case Ui is 
1 + ^1+^2 + ^3 an( l this expression shows that Ui does not produce IR singularity (Ui — > 1 
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for ti — > 0). With a simple calculation, we obtain from Eq. 



t ba =h, t bl =t 2 t 3 , t b2 =ht 3 , 

Z* = {b = (1, 0, 0), 6i = (0, 1, 1), b 2 = (1, 0, 1)}, 

A bo = {xi(l, 1,0) + x 2 (l, 0,1) + x 3 (0, 1,0) + 2:4(0, 0,1) 11, a: 2 , 13, 2:4 >0}, (8) 
A 6l = { Xl (l, 0,0) +3:2(1, 1,0) +s 3 (l, 0,1) I ari,X2,x 3 > 0}, 
A b2 = 0. 

The sub-domains are shown in Fig. [2] The fact that the last sub-domain is empty means 
that term £1^3 never becomes dominant. The second domain is a triangular cone for which 
variables x take simple range of values. However, the first sub-domain includes four param- 
eters in 3-dimensional space. When this sub-domain is divided further into two triangular 
cones, this redundancy disappears: 

A«UAf\ 

On Oo ' 

{a:i(l,l,0) + 12(1, 0,1) + 2:4(0, 0,1) I x u x 2 ,x 4 > 0}, 
{a?i(l,l,0) + as s (0,l,0)+x 4 (0,0,l) | x u x 3 ,x 4 > 0}. 

Now we subsequently change variable from yj to Xi and then to Zi = e~ Xi . The integration 
domain for z becomes 3-dimensional unit cube. Finally we obtain a sector decomposition 
as the following: 

sector variables Z^ 1 

A?^ h = Z!Z 2 , t 2 = z 1: t 3 = z 2 z 4 z 1 z 2 (si2 - s 23 z 4 - s 4 z 2 z 4 ) 
Al o > ti = zi, t 2 = Z1Z3, t s = z A z-i{-sx2 - s 23 z 3 z 4 - S4Z4) 

A bl ti = ZiZ 2 Z 3 , t 2 = Z 2 , t 3 = Z 3 Z 2 Z 3 (-Si 2 Zi - 5 2 3 - S 4 Z!Z 3 ) 

It is easy to calculate Jacobian, which is found to be a monomial of z. 

For the general case, the decomposition of the integration domain in accordance with 
the classification of the asymptotic behavior, corresponding to Eqs. (J6j) and ([8]). becomes a 
problem in Euclidean geometry. It is not so hard to solve this when one uses convex hull al- 
gorithms developed in computational geometry. This decomposition is uniquely determined 
once a polynomial is given. It is independent of the choice of an algorithm. In order to 
obtain final representation of sector decomposition, corresponding to the previous table, it 
is necessary to triangulate convex polyhedral cones. This is performed with triangulation 
algorithms. The triangulation of convex polyhedral cones is not determined uniquely and 
the number of sectors will depend on the algorithm. 

4 Test implementation and conclusion 

We have made a test implementation of this method. Our convex hull algorithm is a modified 
one for convex polyhedral cones based on the algorithm described in Ref. [TT] for polytopes. 
For triangulation of convex polyhedral cones, we have developed our own algorithm. The 
input to our program is given by the package described in Ref. |12j and output is passed 
to the same package for the separation of poles and integrations of coefficients. We have 
checked our program by comparison with another package qhull (Ref. |13| ) for the convex hull 
algorithm, calculating integration volume for the triangulation, and comparison of integrated 
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Table 1: Comparison of number of sectors among different methods. Numbers in column 
"H" are cited from Ref. [I]. Columns "A", "B", "C", "S" and "X" indicate corresponding 
strategy described in [5] and [7] . As shown in [7, , "F" means that the sector decomposition 
fails and "M" means that the memory overflow happened during the sector decomposition 
on a 8Gb machine. The numbers with "*" are given by [14]. "This method" indicates the 
number of sectors obtained by our method. "Exponential S.D." indicates the number of 
sectors before the triangulation. 
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values of several diagrams with references. The number of decomposed sectors are shown 
by Table [U 

We have proposed a factorization algorithm in sector decomposition. Our method is 
based on a classification of asymptotic behavior of polynomials. The problem is converted 
to one for convex bodies in Euclidean space. In order to find sector decomposition for a given 
integrand, we employed algorithms developed in computational geometry. This method is 
deterministic and never falls into an infinite loop. A test implementation shows that the 
number of decomposed sectors is less than iterated sector decomposition combined with 
several strategies. 

This work is supported in part by Ministry of Education, Science, and Culture, Japan under 
Grant-in- Aid Nos. 20340063 and 21540286. 
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